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ABSTRACT 

Context. In recent years, day-side emission from about a dozen hot Jupiters has been detected through ground-based secondary eclipse 
observations in the near-infrared. These near-infrared observations are vital for determining the energy budgets of hot Jupiters, since 
they probe the planet's spectral energy distribution near its peak. 

Aims. The aim of this work is to measure the Kj-band secondary eclipse depth of WASP-33b, the first planet discovered to transit an 
A-type star. This planet receives the highest level of irradiation of all transiting planets discovered to date. Furthermore, its host-star 
shows pulsations and is classified as a low-amplitude S Scuti. 

Methods. As part of our GROUnd-based Secondary Eclipse (GROUSE) project we have obtained observations of two separate 
secondary eclipses of WASP-33b in the K,-band using the LIRIS instrument on the William Herschel Telescope (WHT). The telescope 
was significantly defocused to avoid saturation of the detector for this bright star (K~7.5). To increase the stability and the cadence 
of the observations, they were performed in staring mode. We collected a total of 5100 and 6900 frames for the first and the second 
night respectively, both with an average cadence of 3.3 seconds. 

Results. On the second night the eclipse is detected at the 12 -cr level, with a measured eclipse depth of 0.244;^55q5o This eclipse 
depth corresponds to a brightness temperature of 3270^1^^ K. The measured brightness temperature on the second night is consistent 
with the expected equilibrium temperature for a planet with a very low albedo and a rapid re-radiation of the absorbed stellar light. For 
the other night the short out-of-eclipse baseline prevents good corrections for the stellar pulsations and systematic effects, which makes 
this dataset unreliable for eclipse depth measurements. This demonstrates the need of getting a sufficient out-of-eclipse baseline. 

Key words, techniques: photometric ~ stars: individual: WASP-33 ~ planets and satellites: atmospheres 



1. Introduction 

In recent years, there have been many measurements of thermal 
, emission from the atmospheres of hot Jupiters, especially in the 
i mi d-infrare d using the Spitzer Space Telescope (e.g. the review 
■ by|D eming|l2009>) . These Spitzer observations probe the thermal 
, emission of hot Jupiters redward of the peak of their spectral 
' energy distribution (SED), and thus measure the planet's light in 
' the Rayleigh- Jeans tail of their emission spectrum. 

Observations in the near-infrared, on the other hand, 
typically probe the planet's emission spectrum around or 
even blue-ward of its peak, and therefore provide essential 



information on the planet's total energy budget. During the past 
three years several measurements of planetary emission short- 
ward o f 2.5A(m have been obtained (de Mooii & Snellea 
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Rogers etal 2009; Anderson et al. 2010: Alonso et aL .2010; 


Gibson etal J 12010: Croll et all 1201 Oal: Lonez-M 
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201C 


;ICroll et al.ll2010h,201 iVde Mooii et al.ll2011 : 


Smith et al. 


2011 


; Caceres et al.l 2011: iDemine et al. 2012f). most of these 



are in the K-band (A— 2.15pm) where the planet-to-star con- 



* Photometric timeseries are only available in electronic form at the 
CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat ?J/A+A/| 



trast is most-favourable for observations with ground-based 
telescopes through the available atmospheric windows. 

From the combination of the measurements at multiple 
wavelengths, a picture is emerging that there are (at least) two 
types of hot Jupiter atmospheres, those which show a thermal 
inversion, and those which do not. It has been proposed that the 
presence of the inversion layer is set by the stellar irradiation, 
where at high levels of irradiation the planet's stratosphere is 
hot enough to keep a strongly absorbing compound in the gas 
phase, while at lower irradiation levels the com pound condenses 
out and disappears from the gas-phase (e.g. Burrows et al.ll2007l: 
iFortnev et al.ll2008l) . 

iKnutson et alJ ( 1201 C f) proposed an alternative scenario for 
the presence or absence of a strong absorber in the highest layers 
of the planetary atmosphere. In their scenario the absorber can 
be destroyed by strong UV emission from the planet's host-star, 
due to stellar activity. For higher levels of stellar activity, which 
result in a higher UV flux, the absorbing compound is possibly 
removed, resulting in a non-inverted atmosphere. Note that the 
inference of an inversion layer has recently been questioned by 
Madhusu dhan & Seagej (l2010l) . who point out that for several 
planets there is a degeneracy between the atmospheric temper- 
ature structure and the chemical composition of the planet's at- 
mosphere. 



1 



De Mooij et al.: The GROUSE project III: The secondary ecUpse of WASP-33b 



In this paper we present observations of two secondary 
eclipses of the very hot Jupiter WASP-33b in Kj-band. These 
are part of the GROUnd-based Secondary Eclipse project 
(GROUSE), which aims to use ground-based telescopes for 
exoplanet secondary eclipse observations in the optical and 
near-infrared. As part of this project we have already pub- 
lished Kj-band detections of the secondary eclipses of TrES-3b 
dlEMooij & Sne llen 2009) and HAT-P- l b (de Mooij et al. 201 1) 

WASP-33b jConier Cameron et al.l 1201 Ol) is currently the 
only known planet to transit an A-type star (Teff=7430+100K) 
which it orbits in ~1.22 days. This makes WASP-33b the most 
irradiated planet known to date, with an irradiation of 1.2-10"' 
erg/sec/cm^. This high level irradiation results in an expected 
day-side equilibrium tempe rature of 3250K. Recent observa- 
tions by ISmith et al] (1201 ll) indeed show a very high bright- 
ness temperature at 0.9//m of 3466±140K. Additionally, since 
the host-star is relatively hot, the expected UV flux it receives is 
also high, making it an ideal candidate to investigate the influ- 
ence of a high UV flux on the temperature structure of a planet's 
atmosphere. 

In addition to being the first transiting planet discovered to 
orbit an A-type star, WASP-33b is also t he first planet to transit 
a puls ating star. In the discovery paper, I CoUier Cameron et alJ 
(l20Toh find evidence for non-radial pulsations in their spectral 
time-series, and tentatively classified WASP-33 as a y Doradus 
pulsator, which is a class of non -radial pulsators with period s 
of ~0.3 d ays or longer (see e .g. [Handler & Shobbroo kl l2002l) . 
Recently, iHerrero et al.l (1201 ll) analysed photometric time se- 
ries for WASP-33 and found a pulsation period of 68.6 min- 
utes, which, when converted to the pulsation parameter Q, the 
product of the pulsation period an d the square-root of the mean 
stellar density (e.g. [Breger 199CS iHandler & Shobb rook 2002), 
is comparable to that of 6 Scuti stars, and well outside the range 
of y Doradus stars. The observed stellar pulsations have a mea- 
surable impact on the transit a nd eclipse m easurements for this 
plane t (e.g. .Herrero et al.ii20"Tl1: ISmith et ai.j.201 k.Deming et alj 
I2012h . 

In Sect. |2]we present our observations and data reduction. In 
Sect. [3]stellar pulsations and the light cui-ve fitting are presented. 
Subsequently we discuss the results in Sect. |4l and finally we 
will give our conclusions in Sect. |5] 

2. Observations and data reduction 

2.1. Observations 

The secondary eclipse of WASP-33b was observed on two 
nights, on August 18, 2010 and September 20, 2010, in the 
Kj-band with the Lo ng-slit Intermediate Reso lution Infrared 
Spectrograph (LIRIS; lAcosta-Pulido et aL I I2OOI instrument on 
the William Herschel Telescope (WHT) on La Palma. 

The pixel scale of LIRIS is 0.25 arcsec per pixel, yielding 
a field-of-view of 4.2 by 4.2 arcminutes, large enough to ob- 
serve both WASP-33 and a reference star of similar brightness 
simultaneously. Since WASP-33 is very blight, exposure times 
of 1.5 seconds were used in order to avoid saturation of the de- 
tector. As an additional measure to prevent saturation, the tele- 
scope was strongly defocused. This is a well proven strategy 
also used for other GROUS E observations (Ide Mooii & SnellenI 
120091: Ide Mooii etalj|2011h . which should also reduce the im- 
pact of flat-field inaccuracies by spreading the light over many 
pixels, thereby minimizing the impact from uncorrected pixel- 
to-pixel sensitivity variations. To keep the observations as stable 
as possible, and in order to reduce the cycle time, the observa- 
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Fig. 1. Raw lighcurves for WASP-33 and the reference star (mul- 
tiplied by 2.08 for plotting purposes) for the night of August 18, 
2010 (top panel) and for the night of September 20, 2010 (bot- 
tom panel). The vertical dashed lines indicate the expected be- 
ginning and end of the targeted eclipse. The solid, grey, lines 
show the airmass during the nights, scaled to match the stellar 
flux in the first hour of the observations. 



tions were performed in staring mode. Since this method does 
not allow for background subtraction using the science frames, 
a set of sky frames were obtained after the observations on both 
nights for sky- subtraction purposes. 

On August 18 (Night I) the observations started at 00:45 UT 
and lasted for ~4.5 hours. The weather conditions during the 
night were photometric, as can be seen from the raw light curves 
shown in the top panel of Fig. [Tj A total of 5 100 science frames 
were obtained with an average cadence of 3.3 seconds. The first 
three frames of a sequence of framed are known to suffer from 
the reset anomaly, which is seen as an anomalous structure in the 
background. These frames are therefore excluded from further 
analysis, which results in a total of 4994 frames. 

The observations on September 20 (hereafter night II) were 
taken between 22:35 UT and 05:00UT. During the first part of 
the observations the conditions were photometric, however dur- 
ing the last few hours occasional clouds moved across the image, 
absorbing up to 65% of the light (see Fig.[TJ. A total of 6693 sci- 
ence frames were obtained with an average cadence of 3.3 sec- 
onds per frame, excluding 207 frames due to the reset-anomaly. 



' The observations consisted of multiple sequences of 100-200 
frames 



2 



De Mooij et al.: The GROUSE project III: The secondary ecUpse of WASP-33b 



2.2. Data reduction 

The data-reduction for both nights was performed in the same 
way. All frames were corrected for crosstalk along rows of the 
detector, which is present at a level of 10"^ of the total flux along 
the rows of all four quadrants. Subsequently we performed a 
non-linearity correction on all the frames using our own non- 
linearity measurements which were created from a set of dome- 
flats at a constant level of illumination but with varying exposure 
times. After these corrections the images were flat-fielded using 
a flat-field created from bright and dark twilight flats. 

A background map was constructed from the set of dithered 
images obtained after the eclipse observations. These images 
were reduced in the same way as the science images, and, af- 
ter filtering out the discrepant pixels in time to remove (faint) 
stars, were subsequently combined. The resultant background 
map was then scaled and subtracted from the individual science 
images. 

After background subtraction, aperture photometry was per- 
formed on both WASP-33 and the reference star using an aper- 
ture of 1 8 and 26 pixels for night I and night II respectively. Any 
residual sky background was determined in annuli between 30 
and 50 pixels for night I and between 40 and 60 pixels for night 
II. The flux in the annuli was clipped at 5cr to avoid outliers in the 
background (such as hot pixels) from affecting the data. Finally, 
the light curve of WASP-33 was normalised with that of the ref- 
erence star, and the resultant light curves for the two nights are 
shown in the top panels of Figs.|3]and|4] 



3. Correction for systematic effects and stellar 
pulsations 

3. 1 . Stellar pulsations 

IColher Cameron et al.1 (12010') noted th at the host-star o f 
WASP-33b is a pulsator. Observations bv iHerrero et al.l ( 201 1 



indicated a dominant period of 68.5 minutes. iSmith et al.l (1201 1 



observed a secondary eclipse of WASP-33b in a narrowband fil- 
ter at 0.91;um, and found three pulsation periods in their data, 
at 53.62, 76.52 and 41.85 minutes, all wit h amplitude between 
0.4 mmag and 0.9 mmag. Observations bv lDeming et al.l (1201 2h 
showed various pulsation periods, 68 minutes for observations of 
two separate events in the mid-infrared from the Spitzer Space 
Telescope, 7 1 minutes for one ground-based event in the K-band, 
146 minutes for observations in the J-band and 54 & 126 min- 
utes for a second set of K-band observations. For all their periods 
the amplitudes were larger than 1 mmag, with the two nights of 
K-band data showing amplitudes in excess of 2 mmag. 

For our data, the stellar pulsations are clearly visible in the 
light curve for night I, while for night 2 the variability is less 
apparent. In order to determine the period(s) of the stellar pulsa- 
tions, a periodogram of the light curves was created. Since the 
planetary eclipse signal and possible systematic effects can influ- 
ence the periods found in the data, both a scaled eclipse model 
as well as a model for the systematics based on instrumental ef- 
fects (see section [J!2T l were fitted to the data using a simple linear 
regression algorithm and subsequently divided out before deter- 
mining th e periodogram. The generahsed Lomb-Scargle formal- 
ism from IZechmeister & Kiirsterl (l2009h was used to construct 
the periodogram, and the results for both nights are shown in 
Fig.|2] 

For night I the strongest peak in the periodogram is found at 
a period of 64.5 minutes, and there is a weaker peak at 43.2 min- 
utes. For night II, there are four peaks visible, the strongest peak 




1.0 

Period (hours) 

Fig. 2. Normalised periodograms of the light curves for the two 
separate nights. The top panel shows the periodogram for night 
I and the bottom panel shows the periodogram for night II. The 
dashed lines in both panels indicate the periods used in the fit- 
ting. 



is found at 52.1 minutes, while three weaker peaks are found 
at 43.3, 65.3 and 83.9 minutes. The signals at all the periods 
in both datasets have a false-alarm probability, estimated using 
Monte-Carlo simulations, of below 0. 1 %, indicating that the pe- 
riodicities are very likely real, although they might not be as- 
trophysical in origin. In both datasets we find a periodic signal 
at ~65 m inutes, which difi^ers from the period of 68.5 minutes 
found by iHerrero et al.l (1201 \\ and the 68 minute period from 
the two Spitzer datasets from Deming et al. (2012). However, the 
short time span covered during each night is not sufficient to get 
a very tight constraint on the period, and therefore the periods 
could be consistent with 68 minutes. The ~43 minute period is 
seen in the measurements from both our nights, as well as in 
the data from Smith et al.l (1201 ih . The period around ~52 min- 
utes is fo und in both our sec ond night of data as we ll as in the 
data from lSmith et all (1201 ll) and lDeming et al.1 (1201 2 ). although 
there it is not the dominant frequency. We caution, however, that 
the periodograms used for the frequency analysis were created 
with data that was only partially corrected for systematic effects, 
and therefore can still be influenced by residual (quasi) periodic 
systematic effects. 

3.2. Light curve fitting 

Since the Ught curve is the result of a combination of three ef- 
fects, the stellar variability, systematic effects related to both the 
instrument and the Earth's atmosphere and the secondary eclipse 
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Table 1. Fitted parameters and uncertainties from the MCMC analysis of the hghtcurves of WASP-33b. 



Parameter 


Night 1 




Night 2 






unit 


Instrumental parameters 


Polynomial 


Instrumental parameters 


Polynomial 


F„/F. 


(%) 


0.140+0.007 


0.092+0.017 


0.245+0.009 


0.246+0.018 


X 




0.06±0.01 




-0.12+0.02 




airmass 




-0.76±0.03 




0.34+0.03 




sky 




0.27±0.03 




-0.53+0.03 




Cl 






0.42+0.01 




0.10+0.01 


C2 






-0.09+0.02 




-0.01+0.02 


C3 






0.10+0.01 




-0.13+0.01 


Pi 


(minutes) 


63.95±0.39 


62.00+0.51 


52.65+0.39 


52.28+0.51 


P2 


(minutes) 


43.22 (fixed) 




43.34 (fixed) 




P3 


(minutes) 






65.29 (fixed) 




P4 


(minutes) 






o3.9U (tixea) 




Al 


(%) 


0.095±0.004 


0.082+0.005 


0.056+0.006 


0.058+0.006 


A2 


(%) 


0.041+0.004 


0.042+0.004 


0.017+0.005 


0.014+0.005 


A3 


(%) 






0.011+0.006 


0.011+0.006 


A4 


(%) 






0.031+0.006 


0.013+0.006 



Notes. Baseline coefficents for the instrumental model are the x-position, airmass and sky level and for the polynomial baseline the coefficients 
are Ci to C3. Pi to P4 are the periods used for the stellar pulsations, and Ai to A4 are the corresponding amplitudes. 



of WASP-33b, a fit for all three effects is performed simultane- 
ously. 

For the stellar pulsations the period of the dominant mode 
is left as a free parameter, although with a penalty for the 
of the form (P-Po)^/cr^, with Pq the period determined from the 
periodogram, and o-p set to 1 minute. The periods of the other 
modes was kept fixed to periods found in the analysis of the pe- 
riodogram as described in the previous section. For all the modes 
the offset in phase and the amplitude of the pulsations were al- 
lowed to vary freely. 

For the fitting of systematic effects two diff'erent methods 
were used. For the first method, the systematic effects are con- 
sidered to be due to the change of position on the detector, 
the airmass and the difference in sky background between the 
two quadrants. This is similar to wh at was used in the previous 
papers from t he GROUSE project (de Mooii & Sn ellen"2009^. 
Ide Mooii et al., 20 11) . In the second method, the systematic ef- 
fects are mode l led usi ng low order polynomials, as also used in 
Ide Mooii et alj (1201 2h for the near-infrared transit observations 
ofGJ1214b. 

The secon d ary eclipse was modelled using the 
[Mandel & Asol' ("2002) formalism. We used the parameters 
from Collier Cameron et al. (2010) for the impact parameter, 
semi-major axis, orbital period and planet-to-star size ratio, 
while the orbit of the planet is assumed to be circular. This 
assumption is reasonable since the planet orbits extremely close 
to its host-star, which should result into a rapid damping of 
eccentricity. In addition, on the second night an eclipse-shaped 
dip in the light curve centered on (p ~Q.5 is readily visible (see 
the right panel of Fig. [3]l. 



Before fitting, outliers were removed by excluding all points 
that were more than 0.9% away from a median smoothed light 
curve with a box size of 5 1 points, as well as all points for which 
the flux of the individual stars, corrected for airmass, dropped 
below 90%. In this way a total 21 and 414 points were excluded 
during the first and second night respectively. In addition, there 
is a feature present in both light curves at the same time after 
the start of the observations (after 0.1795+0.0025 days) that is 
not at an identical point during the planet's orbit and therefore 
most likely due to an, as yet, unidentified instrumental effect. 
Excluding all the points that were obtained during this feature 
removes an additional 128 frames. 

The light curves were fitted with 9 free parameters (1 for the 
eclipse, 3 for the systematic eff'ects and 5 for the stellar pulsa- 
tions) for night I, and 13 free parameters for night II (due to 2 
additional periods found in the data). The two nights were fit- 
ted separately using a Markov-Chain Monte Carlo method. Per 
night, 5 sequences of 2 million steps were generated, trimming 
the first 200,000 points to avoid any contamination from the ini- 
tial conditions. The c hains were combined after checking that 
they were well mixed dGelman & Rubin. 1992) . 



4. Results and discussion 

The best fit values of the eclipse depth and their fornal uncer- 
tainties for night I are 0. 140+0.007% and 0.092+0.017% for the 
fit with instrumental parameters and polynomials respectively, 
while for night II the best-fit eclipse depths are 0.245+0.009% 
and 0.245+0.018% for the two respective anlyses. These results 
are given in Table [1] 
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Fig. 3. Light curves for the secondary ecHpse of WASP-33b for the night I (left panels) and night II (right panels). Top panels: the 
light curves of WASP-33 normalised with those of the reference star, overplotted is the best fitting 'full' model with a low order 
polynomial baseline coiTection, stellar pulsations and the eclipse. Middle panels: The light curves corrected for the trends in the 
baseline and stellar pulsations, clearly showing the transit. Bottom panels: The residuals after subtracting the best-fit model. The 
thick points with errorbars in these figures show the data binned by 50 points. The vertical dashed lines show the expected times for 
first to fourth contact. 



The differences found between methods in the first night are 
significantly larger than the uncertainties in the eclipse depth as 
estimated from the MCMC analysis. First of all the first night 
suffers from a strong peak, possibly due to stellar pulsations, 
right in the middle of the eclipse. We attribute this to a prob- 
lem with the observations on night I. The relatively short out-of- 
eclipse baselines available for the first eclipse observation ham- 
pers the removal of the systematic effects from the stellar pulsa- 
tions, as well as from the instrumental and atmospheric effects. 
This is clearly illustrated when looking at the correlations be- 
tween the parameters used for the removal of the systematics 
and the eclipse depth, as shown in Figs. IA.11 to lA.4l We there- 
fore conclude that the first night of data is no usable for a reliable 
eclipse measurement. 

To assess the impact of correlated noise, we redid the anal- 
ysis for night II after binning the data by 50 points (~3 min- 
utes). Although overall the parameters are the same, we find 
larger uncertainties in the eclipse depths with the best-fit values 



of 0.255+0.028% and 0.242+0.035% for the fit with insti'umen- 
tal parameters and polynomials respectively. To assess the im- 
pact of (uncorrected) red noise on the measured eclipse depths 
in another way, th e residual permutation method was used (e.g. 
iGillon et alJ2"007h . The best fit model is subtracted from the light 
curve, and these residuals are then shifted by n points, wrap- 
ping the light curve around, so that the points that are shifted 
beyond the end of the lightcurve are inserted at the beginning. 
The best fit model is then added back to the data, and this new 
light curve is fitted again. The interval between 16% and 84% 
of the distribution of the best-fitting eclipse depths is used for 
the 1-cr uncertainties on the eclipse depth. To speed up the resid- 
ual permutation analysis, instead of adding back the full model, 
which includes the stellar pulsation, trends in the baseline and 
eclipse depth, we only used the trends in the baseline and eclipse 
depth, since the correlation between the parameters for the stel- 
lar pulsations and the eclipse depth is weak. From the residual 
permutation analysis we also find larger uncertainties for both 
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Fig. 4. Same as Fig.|3]but the systematic effects are now modelled by Linear relations with instrumental parameters. 



decorrelation methods, with eclipse depths of 0.244;^q Jjjo'^ for 
a baseline fitted with instrumental parameters and 0.249 ^Q g^j*^ 
for a polynomial baseline fit. In all cases the uncertainties are 
higher than for the MCMC analysis of the unbinned data but 
comparable to the MCMC analysis of the binned data, which is 
expected in the presence of red noise. 

Since there is a strong correlation between the coefficients 
for the polynomial baseline fit and the eclipse depth (see 
Fig. IA.2b . we use the fit of the baseline with instrumental pa- 
rameters for the remainder of the paper, since the correlation 
between different parameters is much weaker. We note that the 
polynomial baseline correction for this night gives the same 
eclipse depth, however with a larger uncertainty. 

The measured eclipse depth of 0.244^^ ^20 corresponds 
to a brightness temperature in the K^-band of 3270;^ [g^ K. 
This brightness temperature was calculated using t he solar- 
metallicity NextGen models dHauschildt et alj [l99^ interpo- 
lated to the stellar parame ters of WASP-33 determined by 
ICoUier Cameron et alJ (1201 Oh (Teff=7430K, Iog(g)=4.294). 

Currently there are several other measurments of the sec- 
ondary eclipse of WASP-33b, Smith et al. (2011) obstained 
data in a narr owband filter at ~9 100A, showing a depth of 
0.09+0.016%. iDeming et aP (l2012h observed both in the Ks- 



band and in the IRAC 3.6jum and 4.5fim bands, obtaining 
eclipse depths of 0.27±0.04%, 0.26+0.05% and 0.41+0.02% 
respecively. These brightness measurements correspond to 
a brightness temperature of 3490+140 K, 3415+130 K, 
2740+225 K, 3290+100 K for the SII, K„ IRAC 3.6jum and 
4.5fj.m respectively. Although a planet's brightness temperature 
can be a strong function of wavelenght, most of the measure- 
ments are consistent at the Icr level. 

As can be seen in the middle right panel of Fig. ID the eclipse 
appears to end earlier than expected from the model. Although 
systematic effects are the most likely cause of this, it is worth 
noting that a narrower width of the eclipse would indicate that 
the orbit of WASP-33b is eccentric. If this is the case, by com- 
bining the ratio between the transit and secondary eclipse du- 
rations with the time of mid-eclipse, a direct measurement of 
both the e ccentricit y and the argument of periastron is possi- 
ble (e.g. lCharbonneau et aLll2005 ). Although a full fit is beyond 
the scope of this work, we can estimate the change in duration 
from the light curve. The duration of the secondary eclipse is 
shorter than the transit duration by 0.01 in phase, such that the 
eclipse duration corresponds to ~90% of transit duration. From 
this we estimate esin(w)~0.05. Since the ingress appears to be 
at the expected time, the time of mid-eclipse is in this case 
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Fig. 5. Spectral energy distribution of WASP-33b. Top panel: 
Ecli pse depths in K^-ba nd (this work) and in the Sllo.gi^m-filter 



from lSmith et al.l(l201 1 ). Bottom panel: brightness temperatures 
in the two bands. Overplotted in both panels are the expected 
eclipse depths/brightness temperatures for a zero-albedo homo- 
geneous day-side (solid line), for an instantly re-radiating day- 
side (dashed line) and for the best-fit effective temperature (dot- 
ted line). 



also slightly earlier than expected. From the shift we estimate 
ecos((y)~0.008. Combining these two estimates we find an ec- 
centricity of ~0.05. We again caution that systematics can easily 
give rise to an apparent non-zero determination of the eccentric- 
ity, which is for instance seen for the secondar y eclip se of TrES- 
3b (Ide Mooii & Snellen! 120091: iFressin et"an l20Tot ICroll et al.1 



I2010bh . We therefore do not advocate this non-zero eccentric 
ity scenario based on these data. 

4.1. A low albedo and rapid re-radiation of incident light 

With the exception of the Spitzer 3.6;um measurement from 
[Oeming et al. (2012), all of the currently available eclipse mea- 
surements for WASP-33b point towards a very hot day-side tem- 
perature. If we assume that the measured brightness tempera- 
tures are representative of WASP-33b's equilibrium temperature, 
and are not generated deep inside the planets atmosphere, where 
the temperatures are even higher, we can constrain the planet's 
equilibrium temperature to Te//p=3298^g^K (see also Fig. [S]). 
Note that the basic assumption that the brightness temperature 
equals the eff'ective temperature does not necessarily has to be 
the case, but detailed modelling of the available measurements 
is beyond the scope of this paper. From this equilibrium temper- 
ature we can further constrain the re-radiation factor and albedo. 
For the re-radiation factor we used the f description, as used by 




Fig. 6. Equilibrium temperature of WASP-33b for different 
albedo (A) and re-radiation factors (f). Solid contours show 
lines of constant temperatures at 150 K intervals. Overplotted 
are lines of constant temperature for the measured brightness 
temperatures (dashed lines), labeled with the bandpass they 
were obs erved (D12 indicates the K^-band measurement from 
Deming e t al.l ( l20I2h . The line labeled "Avg." indicates the line 
for constant temperature of the effective temperature determined 
from the literature measurements. Vertical (dashed) lines indi- 
cate the re-radiation factors for a homogeneous day-side temper- 
ature (f=l/2), and for an instantly re-radiating day-side (f=2/3). 



iLo pez-Morales & SeaeeddloOTl ): 

/R \'/^ 



(1) 



With R» the stellar radius, a the semi-major axis, Ag the bond 
albedo and 1/4 < f < 2/3, where f=l/4 is for a homogenous tem- 
perature distriubution across the planet and f=2/3 is for instant 
re-radiation. 

In Fig.|6]we show a simple model of the equilibrium temper- 
ature as a function of albedo and re-radiation factor. In addition 
we show lines of constant effective dayside temperature with the 
observed brightness temperatures and derived effective tempera- 
ture overplotted. As can be seen the measurements require a very 
low albedo and a very short re-radiation time scale, such that all 
the stellar flux is absorbed and rapidly re-radiated without hav- 
ing time to advect to the night-sid e of the planet. This is consis- 
tent with the findings of ICowan & Agol (201 1), who study the 
albedo and redistribution efficiencies for a large sample of hot 
Jupiters, and find that the hottest planets (in their sample) have 
a low albedo and a low efficiency of the advection of absorbed 
stellar flux to the planet's night side. 

The low redistribution efficiency suggests that the re- 
radiation timesc ales are short, and tha t the planet prob ably has an 
inversion layer (iFortnev et al.ll2008l) . iKnutson et al.l ((2010) hy- 
pothesise that an increase in the UV-flux from an active star can 
cause a shift in the photochemistry such that the efficient ab- 
sorber is removed from the gas-phase. The high incident UV-flux 
on WASP-33b would argue against this. It should be noted that 
for active stars most of the UV flux is emitted in the Lyman a 
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line, while for WASP-33 it is likely that the UV continuum emis- 
sion dominates. To investigate the influence of the UV-radiation, 
photo chemical modelling will be necessary (e.g. iZahnle et ah 
120091) . 

5. Conclusion 

We have presented our results of K^-band observations of the 
secondary eclipse of WASP-33b, the most irradiated planet 
known to date. The measured eclipse depth for the second night 
is 0.244;^g g2g %, which results in a brightness temperature of 
3270;^ [gg K. This high brightness temperature, if representative 
for the planet's equilibrium temperature, requires a very low 
albedo and a high (f>0.5) re-radiation factor 

Combining our K, -band measurement with the measurement 
of lSmifli etaiT(l20TTl) . we can fit a simple blackbody function 
to the spectral energy distribution, and determine an effective 
temperature of Teff,p=3290+^^K. 

We also find that stellar pulsations of the 6 Scuti host-star, 
WASP-33, appears to have switched modes between the two 
nights, which are located a month apart, and also differ from 
the measurements by Herrero et al. (2011). We caution, how- 
ever, that this could be due to systematic effects which could 
also have strong periodicities. 

The measurements on the first night suffer from strong resid- 
ual systematics and stellar pulsations that cannot be fully cor- 
rected due to the short out-of-eclipse baseline, and demonstrates 
the need for observing a target for as long as possible. 
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Fig. A.l. Correlations between different parameters from the MCMC analysis for night I. The baseline is modelled using a polyno- 
mial with coefficients ci to C3. The correlation between the parameters and the ampUtudes of the stellar pulsation at lower signifi- 
cance are not shown, since their impact on any of the measured parameters will be minimal. The amplitude for stellar pulsation is 
for a pulsation period of 62.00 minutes. 
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Fig. A.2. same as Fig. lA.ll but now for night II. The amplitude for stellar pulsation is for a pulsation period of 52.28 minutes. 
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Fig. A.3. Correlations between different parameters from the MCMC analysis for night I. The baseline is modelled using the ob- 
served instrumental parameters. The correlation between the parameters and the amplitudes of the stellar pulsation at lower signif- 
icance are not shown, since their impact on any of the measured parameters wiU be minimal. The amplitude for stellar pulsation is 
for a pulsation period of 63.95 minutes. 
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Fig. A.4. same as Fig. lA.3l but for night II. The amplitude for stellar pulsation is for a pulsation period of 52.65 minutes. 
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